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Fluctuations of entropy production in driven glasses 
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We study by Montecarlo simulation the non-equilibrium stationary behavior of a three-dimensional 
stochastic lattice gas with reversible kinetic constraints and in diffusive contact with two particle 
reservoirs at different chemical potential. When one of the boundaries is placed in the "glassy 
phase", the statistics of entropy production fluctuations driven by the chemical potential gradient 
satisfy the Gallavotti- Cohen fluctuation theorem provided that an "effective" entropy production is 
introduced. 
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Introduction. The Gallavotti-Cohen fluctuation the- 
orem (FT) is a result of remarkable generality concern- 
ing the distribution of entropy production fluctuations 
over long time intervals in stationary open systems. It 
states that the ratio of the probabilites of observing a 
given average entropy production er r , over a time inter- 
val t, to that of observing the opposite value — ay is 
e TOV Q. It has been recently shown Q that this result 
implies the fluctuation-dissipation theorem (FDT) and 
Onsager's reciprocity relations in the limit of vanishing 
driving forces (i.e., in equilibrium). Therefore, the FT 
appears as an important stepping stone towards the gen- 
eralization of thermodynamics to situations out of equi- 
librium. 

The FT was inspired by a pioneering numerical ex- 
periment on a thermostated fluid under shear flow ||. 
It has been further confirmed by the numerical simula- 
tion of an electrical conduction model |Q, and heat dif- 
fusion in a Fermi-Pasta-Ulam chain [^j, while a turbu- 
lent Rayleigh-Benard convection experiment is still un- 
der analysis Although the theorem was originally 
proved for thermostated Hamiltonian systems driven by 
external forces, under certain "chaoticity" assumptions 
for the dynamics [Q, it also holds for rather "generic" 
stochastic systems endowed with Langevin [JtJ , and diffu- 
sive dynamics ||. 

It is well known, however, that violations of FDT 
appear in a wide class of conservative systems, since 
their relaxational dynamics after a quench in the low- 
temperature or high-density phase is so slow that pre- 
vents them from reaching equilibrium on finite time 
scales ||. In the attempt to formulate a thermology of 
aging systems, the FDT violation has been related to 
the appearance of an effective temperature, T ff = T/X, 
different from that of the thermal bath, T jlO) (for a dif- 
ferent point of view, see [0). The factor X is called 
the fluctuation-dissipation (FD) ratio: it is equal to one 
at high temperatures (where FDT holds), while it may 
become smaller than one for sufficiently long times and 
at low enough temperatures. In view of the connection 



between FDT and FT we can ask ourselves what form 
assumes the latter in those driven systems for which a 
FDT violation is expected in the limit of small entropy 
production. We have therefore considered the following 
question. Given a conservative aging system, C, what is 
the form of FT for the corresponding dissipative system, 
T>, obtained by putting C in contact with two reservoirs 
(one of which in the glassy phase)? The results of our 
computer experiment suggest that for such a driven sys- 
tem a generalized form of FT holds where the entropy 
production is evaluated taking into account the FD ratio, 
X , of the corresponding conservative, undriven, system. 

The driven glass. To corroborate this conjecture we 
consider a three dimensional lattice-gas model of hard- 
core particles kept away from equilibrium by a chemical 
potential gradient (see fig. [l]). At the top and bottom 
layer of the system particles are inserted and removed ac- 
cording the usual Montecarlo rule, mimicking the contact 
of the system with a particle reservoir 1Z± : we randomly 
choose a site on the layer; if it is empty, we add a new 
particle; otherwise we remove the particle with probabil- 
ity e~ /3fl± (/i± > 0). The global effect of the reservoirs is 
to fix the boundary densities at two different values p+ 
and p- , driving a current through the system. A similar 
setup was previously consider by Spohn with the aim 
of studying analytically long-range correlations in a sim- 
ple model of non-equilibrium fluid subject to temperature 
gradients. The link of the microscopic theory with the 
fluctuating hydrodynamics was also investigated, see |lE| | 
and, for a review, |]l6| . Here, we consider a bulk dynamics 
implementing the cage effect in supercooled liquids 
Therefore, the sweeps of creation/destruction of particles 
on the edges are alternated with the following diffusive 
sweep. A particle and one of its neighbouring sites are 
chosen at random; the particle moves if the three fol- 
lowing conditions are all met: 1) the neighbouring site is 
empty; 2) the particle has less than 4 nearest neighbours; 
3) the particle will have less than 4 nearest neighbours 
after it has moved. The rule is symmetric in time and, 
unlike the usual driven diffusive systems where detailed 



1 



balance is broken |p|fl9[ , the system is microscopically 
reversible. 

The dynamical behavior of the system when the two 
reservoirs are at same chemical potential has been investi- 
gated in ]2p|-|22]| . It is worth to recall briefly some results. 
Although the equilibrium thermodynamics of the model 
is trivial, a purely dynamical transition takes place when 
the system is quenched in the glassy phase, namely above 
the threshold value p c ~ 2.0: the density approaches by a 
power-law an asymptotic value p c ~ 0.88, different from 
the equilibrium value allowed by the dynamical evolution 
rule, and the diffusion coefficient of the particles van- 
ishes ro ughly as the inverse of the time elapsed after the 
quench |2(j,|2]J. Consequently the system ages, and the 
mean-square displacement and the conjugated response 
function satisfy a generalized FDT with a FD ratio given 
by X ~ 0.79 §§. 

We are interested here to the case where a boundary 
is in the fluid phase < p c ) and the other one in 
the glassy phase (p + > p c ). The numerical experiment 
is performed in the following way. The system consists 
of a cubic lattice of size L 2 x 2L with L = 10. The 
reservoir 1Z + is located at z — while 1Z- at z = ±L. In 
this way we can take periodic boundary conditions in all 
directions and avoid spurious edge effects p3| . We start 
with a low-density uniform distribution of particles and 
let the system reach a stationary state characterized by a 
time-translation invariant correlation function of current 
fluctuations. We then follow the dynamical trajectory of 
the motion for a time of 10 8 Montecarlo sweeps (MCs) 
along which the observables of interest are evaluated. In 
particular we consider the particle current J(t), defined 
as an extensive quantity (proportional to the transverse 
surface L 2 ) by: 



L-l 



J(t) 



^ J2 b'+(M)-j-(M)], 



(1) 



z=-L 



where j + (z,t) and j_(z,i) are the number of jumps tak- 
ing place at time t through the layer z in the direction 
of, and opposite to, respectively, the externally imposed 
chemical potential gradient. 

Fluctuation theorem and effective entropy production. 
In order to study the statistics of current fluctuations we 
consider the average particle current over a time interval 
of duration r: 



JAt) = - £ J(S), 



(2) 



s=t+l 



which in a steady state does not depend on t. We then 
compute the probability distribution n T {p) of the adi- 
mensional variable 



Jr 



(3) 



where J = lim-r^oo J T . The function tt t (p) also gives 
the distribution of entropy production fluctuations since, 
according the thermodynamics of linear irreversible pro- 
cesses |T^] , the average entropy production a r over a time 
interval r is simply related to the corresponding average 
particle current J T by: 



M-J 



(4) 



Fig. shows the distribution ir T (p) for different values 
of the time interval, r = 1,2,5,10 MCs. The chemical 
potential of the reservoirs is p+ = 2.2 and /i_ = 0. Being 
the system in the "large deviation" regime it can appear 
surprising that the shape of ir T (p), is compatible with a 
gaussian distribution, even though no a priori relation 
is expected with the Sinai limit theorem Jl3| . For a dis- 
cussion see where a similar result was also obtained; 
while a skew distribution was found in 

Once determined the probability distribution tt t {p) of 
entropy-production fluctuations it is simple to check the 
Gallavotti-Cohcn FT, that in our case reads: 



log 



K T (-p) 



-pJ{p+ - p-). 



(5) 



Fig. U shows that the lhs and rhs of the relation (||) are 
linearly related but are not equal. However, an equality 
is obtained by taking into account the "effective chemical 
potential" p e g = Xp (or, equivalently, the FD ratio X) 
of the glassy boundary in the calculation of the entropy 
production: 



OY = J T (p c ff - p-) 



(6) 



When the reservoirs are both placed in the fluid phase, 
where X = 1 (and p c s — p+), the FT should be re- 
covered. The inset of fig. || shows, consistently, that for 
/!+ = 1.8, p- = 0, and r = 1,2,5 MCs the FT is well 
satisfied. 

To rule out that these results are a mere coincidence we 
have also performed extensive numerical simulations with 
different values of the chemical potential In fig. |J, 
we show the results for p_ = -0.4, -0.2, 0.0, 0.2, 0.4, and 
// + = 2.2, with r = 1 MCs, plotted against the effective 
entropy production. We see that also this plot supports 
the conjecture of a generalized form of FT. 

We now give an heuristic argument showing how the 
appearance of the FD ratio here should not be surpris- 
ing. Let us observe that our experimental setup is noth- 
ing else than a device to perform measurements of the 
effective chemical potential, /i c ff, of the system p3. In- 
deed, suppose that we want to measure the value of p c g 
for an aging system in contact with a reservoir 1Z + at 
p+ > p c . We start with the stationary system in contact 
with 1Z+ and 1Z-, with /i_ < p c , the last reservoir acting 
as a "thermometer" . Then, if we disconnect the system 
from 7Z- and let it relax, the density of the correspond- 
ing layer will increase until the average particle current 
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is zero (the value of [i where this condition is first met 
defines /x e ff)- However, since the system is unable to fully 
equilibrate with the reservoir 1Z + (being > /i c ), such 
a density will be lower than p c and therefore the corre- 
sponding chemical potential will be reduced by a factor 
X, i.e., /i c ff = X In principle, this method represent 
a possible way to measure the FD ratio of aging systems 
and gives some insight into the notion of effective tem- 
perature as partial equilibration factor jlC| . 

It is also interesting to observe that the stationary 
density profile can be predicted by a non-linear diffusion 



equation 21 1. However, this matter as well as the validity 
of Einstein relation between diffusivity and conductivity 
will be discussed elsewhere. 

Conclusions and perspectives. We have investigated 
by Montecarlo simulation non-equilibrium stationary be- 
havior of a constrained lattice-gas model driven by a 
chemical potential gradient (driven glass). The ap- 
proach does not involve fictitious thermostating mech- 
anism allowing dissipation to prevent the heating up of 
the system 25 1, neither the question of the interpretation 
of phase space contraction rate as entropy production, 
which has not yet received an unanimous answer (see 
however ^6|). The extension to other systems such as 
driven spin-glasses and to different boundary conditions 
modelling, e.g., a Couette flow, is also feasible pj| . 

In particular, we have studied the distribution of en- 
tropy production fluctuations when one of the reservoirs 
is placed in the glassy phase, /i+ > /i c , and shown that in 
the stationary state a generalized form of FT is satisfied 
which takes into account the effective entropy production 
computed through the FD ratio of the corresponding un- 
driven system. Therefore, the generalized form of FT and 
FDT should be related by the same limiting procedure 
exploited in [^). This vindicates the crucial role played 
by the notion of "effective temperature" in generic, aging 
or stationary, out of equilibrium systems. 

Although the model considered here is much simpler 
than a real supercooled liquid in a temperature gradient, 
it displays interesting phenomena that could be observed, 
hopefully, also in a "true" experiment. 

We conclude by observing that it is not clear what 
happens in a generic driven glass in which the reservoirs 
are both placed in the glassy phase, at different values 
of temperature. In such a strongly non-linear regime of 
small entropy production the correlation of current fluc- 
tuations could exhibit aging and the system never reach a 
stationary state. What would be the form of Gallavotti- 
Cohen FT in such a situation, is presently unknown. 
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FIG. 1. Rayleigh-Benard experiment to study out of equi- 
librium stationary states in driven glasses. The system is cou- 
pled to two particle reservoir at different chemical potential, 
the lower one in the glassy phase, > fi c , and the upper 
one in the fluid phase, /x_ < fi c . Particles flow through the 
bulk from the bottom to the top. 




FIG. 2. Histograms of the probability distribution tv t (p) of 
entropy production fluctuations along a trajectory of motion 
of 10 s MCs sampled at time intervals r = 1, 2, 5, 10 MCs (from 
top to bottom). The chemical potential on the boundaries is 
= 2.2 (glassy phase) and = (fluid phase). The curves 
are all consistent with a gaussian distribution of average 1. 




FIG. 3. Logarithmic probability ratio log(7r T (p)/7r T ( — p)) 
vs. the entropy production over a time interval r, for r = I 
(diamonds), 2 (stars), 5 (squares), 10 (circles) MCs; and 
fi + = 2.2, //_ = 0. The broken line, with slope 1, represents 
the prediction of the FT. The straight line, with slope X, take 
into account the effective chemical potential /u e fi = of 
the glassy boundary. The value of FD ratio, X — 0.79, is 
taken from ref. pj| . Inset: check of FT when both the reser- 
voirs are placed in the fluid phase, fi+ = 1.8 and /x_ = 0, for 
t — 1 (diamonds), 2 (squares), 5 (circles) MCs. 




FIG. 4. Logarithmic probability ratio vs. the effective en- 
tropy production for r = 1 MCs, fx+ = 2.2, and ;U_ = —0.4 
(asterisks), —0.2 (squares), (circles), 0.2 (diamonds), 0.4 
(stars). The straight line with slope 1 represents the predic- 
tion of generalized FT, where fi e g = Xa + and the value of 
FD ratio, X = 0.79, is taken from ref. §21. 
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